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We construct the Euclidean Green functions for the soliton (magnetic monopole) field in the ?7(1)4 
Lattice Gauge Theory with Wilson action. We show that in the strong coupling regime there is 
monopole condensation while in the QED phase the physical Hilbert space splits into orthogonal 
soliton sectors labelled by integer magnetic charge. 
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I. INTRODUCTION 

In this paper we apply the methods introduced in [ 1 ] to the construction of soliton (magnetic monopole) sectors for 
the U(l) Lattice Gauge Theory with Wilson action [2|. The soliton quantization for the U(l) Lattice Gauge Theory 
with Villain action ||, as well as for a large class of models, has been carried out in Jl],[| by constructing the Euclidean 
Green functions of soliton fields as expectation values of suitable disorder operators. These operators are obtained 
by coupling the theory to a generalized external gauge field in an hyper-gauge invariant way, as we shall briefly recall 
below. In the statistical approach this procedure corresponds to the introduction of open-ended line defects. In 
the case of the f/(l) 4 gauge theories such line defects are just magnetic loops carrying a defect (topological) charge. 
Opening up a loop one introduces magnetic monopoles at the endpoints of the current line. An Osterwalder-Schrader 
(O.S.) reconstruction theorem ]|,|| applied to the disorder fields correlation functions permits the identification of 
vacuum expectation values of the soliton field, which can be considered as a charged operator creating magnetic 
monopoles 

We have focused our attention in particular on the vacuum expectation value of the soliton field, which is defined by 
a limiting procedure starting from the two point function: 



Si 



lim 

\x\— >o 



S 2 (0,x) 



The result of our work is that 5*1 is a good disorder parameter for the phase transition occurring in the model [|9|. 
In fact we show that Si is bounded away from zero in the strong coupling regime (/3 <C 1) while is vanishing in the 
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QED phase. The proof of this statement makes use of two different techniques. The strong coupling phase is analysed 
by means of a convergent Mayer expansion, applied to a polymer system obtained by the dual representation of the 
model Besides, the clustering property in the QED phase is proved u sing an adapted version of the expansion in 
renormalized monopole loops originally given by Frohlich and Spencer fll| , |l2|| . The dual representation of the Wilson 
model has a measure given by products of modified Besscl functions: the extimates have been done applying suitable 
bounds on modified Bessel functions for (3 <C 1, while using for (3 large an interpolation to Bessel functions given in 
& 

The paper is organised as follows: in this section we shortly define the disorder fields and correlation functions for the 
U(l) Wilson model and describe their connection with magnetic monopoles. Then we enunciate the reconstruction 
theorem for the U{\) Wilson model. In section || we'll map the model into a polymer system in order to prove that 
Si is nonvanishing in the strong coupling regime and thus the lattice solitons "condense" in the vacuum sector. In 
section [II we'll give the expansion in renormalized monopole loops in order to show that Si = in the weak coupling 



regime: the Hilbert space of the reconstructed Lattice Quantum Field Theory splits into orthogonal sectors labelled 
by the magnetic charge. Finally, in section IV we give some concluding remarks. 



A. Disorder fields and magnetic monopoles 

In what follows we consider all fields as defined on a finite lattice A C Z 4 ; all estimates needed to proof our 
statements are uniform in |A| and therefore extend to the thermodynamic limit (A — > Z ). The partition function for 
the U(l) gauge model is 

z = I vd U fp( d ° P ) . (i) 

J PCA 

where T>0 is the product measure on the 1-form valued in [—it, it), d9 is the field strength defined on the plaquette 
p and (for Wilson action) 

ipp (d9) = e 13 cos{d0) . (2) 

To define a disorder operator we can consider a modified partition function in which an external hyper-gauge field 
strength X is coupled to the dynamical variables: 

Z(X) = f VO Y[ tpp (d9 p + X p ) . (3) 

PCA 

The mean value of the disorder operator is defined by 

and is invariant under the hyper-gauge transformation 

X — > X + d"f, (5) 

with 7 a generic 1-form. The hyper-gauge invariance follows by the redefinition of link variables (6 — * 6 — 7) and 
amounts to say that {V(X)) depends only on the 3-form dX = M . In fact, by Hodge decomposition, on a convex 
lattice one can write 

X = da + 5—M (6) 

and the first term on the right-hand side can be always absorbed in a redefinition of 9. Moreover, M turns out to 
be the dual of the magnetic current density J , because the total field strength G — dd + X satisfies the following 
modified Maxwell equations (in presence of an electric current J E ): 

5G = J E (7) 
dG = M or S *G= J M , (8) 
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where by *G we mean the Hodge dual of the form G. Finally, from (g) follows that the magnetic current is identically 
conserved 

SJ M = . (9) 

Hence we have that a disorder operator as defined by (^) is unavoidably connected to a magnetic current. In this 
language a disorder operator describing an open-ended line defect in the statistical system (jl]) will correspond to 
a current J M describing the birth and the evolution of a magnetic monopole. It is possible to parametrize such a 
conserved magnetic current in a general way as follows: 

J M = 2vr [D - fi] . (10) 

Here f2 is an integer valued 1-form with support on a line (the Dirac string) whose endpoints are the space-time 
locations of monopolcs: 

5Q (x) = H S(x - asi) q, e Z\{0} . (11) 

i 

= (0, -D) is a superposition of Coulomb- like magnetic fields with flux ±qi, spreading out at the time slices where 
the monopoles are located ( ) 

D(x ,x) = ±^-p$^5(x°-x^ (12) 

47T \X — Xi\ 5 

in such a way that SD = SCI. Moreover, hyper-gauge invariance joint to the compactness of the action implies that 
the disorder operator does not depend on the shape of the string. In conclusion we see that a magnetic monopole of 
charge q can be implemented as a defect in the model (|l|) by a 2-form X whose curvature dX plays the role of a dual 
magnetic current. For a monopole- antimonopole pair of charge iq we have 

*J M = dX = 2Trq(*D- *Q) X = 2nq 5 i (*D - , (13) 

where *Q is an integer valued 3-form and *D is the 3-form representing the Coulomb field. 



B. Mixed order-disorder correlation functions 

Performing a Fourier analysis on (^|) we obtain 

z w= e n^M ei(n ' x) ■ ( m ) 

n : Sn— pCA 

rip is the integer valued 2-form labelling the Fourier coefficients. With Ip (n) we indicate the modified Bessel functions 
of order n evaluated in j3 (commonly written as /„(/?)). On the dual form v — *n the constraint Sn = becomes 
dv — and so we can write v = d A and sum over equivalence classes of integer valued 1-forms defined by [A] = {A' : 
dA' = dA}: 

Z{D,tt)=Yl WhidAp) e i2 * q ( A ' D -^ . (15) 

[A] pCA 

Now e l27T <i{ A ,n) = ^ f or integer values of q, and we can conclude that (T>(D,fl)) actually depends only on (xi,qi), 

once we have fixed the shape of the magnetic field D satisfying 5 D — — qiS(x — x{). 

Now one can introduce ordinary fields, preserving hypergauge invariance of expectation values: 

^p(D,fl) = eW> +x 'l . (16) 

The correlation functions to which the reconstruction theorem applies are then given by 

S n ,m {xiqi, ,x n q n ; pi p m ) = (D(xi qi, ,x„q n ) tp Pl ip Pm ) (17) 

for "^2 qt = 0. Correlation functions with non vanishing total charge are defined by the following limiting procedure 
for Qi = Q- 

S n m = lim c q S n+ i, m (xi,qi; ; x n , q n ;x, -q; p x p m ) (18) 

x — *oo 

where c q is a normalization constant. We recall now the reconstruction theorem in the version given in 
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Theorem 1 If the set of correlation functions {<SVi,m} is 

1 ) lattice translation invariant; 

2) O.S. (reflection) positive; 

3) satisfies cluster properties; 

then one can reconstruct from {S„. m } 

a) a separeble Hilbert space TL of physical states; 

b) a vector Q of unit norm, the vacuum; 

c) a selfadjoint transfer matrix with norm \\T\\ < 1 and unitary spatial translation operators fx = 1, ..,d — 1 such 
that 

TQ = UQ = n ; 

d) ft is the unique vector in TL invariant under T and U . 

If moreover the limits ( flf[ j vanish, then TL splits into orthogonal sectors TL q , q £ Z , which are the lattice monopole 
sectors. 

In our case hypotheses 1),2),3) follow from traslation invariance and reflection positivity of the measure denned in 
the standard way from Wilson action [fl4| . In particoular one can easily check the reflection positivity of monomials 
of disorder fields in the dual representation, where they assume a standard form and have support on fixed time slices 
(because the magnetic field spreads out in fixed time planes). 

In the quantum mechanical framework S 2 (x, q; y, —q) is the amplitude for the creation of a monopole of charge q at 
the euclidean point x and its annihilation in y. We are now going to show that in the confining phase the two point 
function S 2 (x, q; y, —q) is uniformely bounded away from zero, while in the QED phase it vanishes at large euclidean 
time distances: 

lim S 2 (x, q ; y, -q) = . (19) 

\x-y\— >oo 

A generalization of these estimates to S n , m implies that in the confining phase there is the so called monopole 
condensation and in the weak coupling region the physical Hilbert space TL decomposes into orthogonal sectors 
labelled by total magnetic charge. 

II. MONOPOLE CONDENSATION IN THE CONFINING PHASE 

In this section we are going to prove monopole condensation in the strong coupling regime: this relies on the fact 
that given arbitrary 7 s]0, 1[, there exists a /3 7 such that for [3 < /3 7 holds 

S 2 (x,q; y,-q) >7 V x,y e A . (20) 

In order to prove our statement we shall adopt the following strategy [[jof : first, we express the disorder field expec- 
tation value in terms of logarithms of the partition functions 

S 2 (x, q ; y, -q) = exp [log Z(D) - log Z] (21) 

and then we prove that [log Z(D) — \ogZ] is close to uniformely in x and y. The main tool we'll use is a cluster 
expansion for [logZ(D)], which we obtain after rearranging Z(D) as the partition function of a polymer gas. 

A. The polymer expansion 

The first step in our program is the polymer expansion for the (modified) partition function of our system (see e.g. 
HQ). Although it is a standard technique we now briefly present its application to the compact U(l) system with 
disorder fields. First note that ZiD, f2) (|f|) can be written summing over closed integer valued 2-forms v as follows: 

z(d,ci) = n a II^M ^ q{A ^ D) , (22) 

v : dv—0 pCA 

with 



4 



1 fi(0) 

In equation (jH) A v is the representative element of the class defined by dA v = v. it is simple to verify that each 
term in the expansion does not depend on the choice of the representative element A v . Moreover we have extracted 
an overall factor rescaling the modified Bessel functions by 1/3(0): the advantage of this choice is that lp(0) = 1 and 
in our expansion around v = we must not carry over tedious factors. 
Let us now give some definitions: the support of a k-form i>w is the following set: 



suppu (fc) = jai € A : x € c k k-ccll with v (k \c k ) ^ oj 



Let us note that following we think of suppv^ as a set of points rather then a set of k-cells. In the same way 
one can define in a natural way the support of a set of k-cells. This permits us to extend to these sets the common 
definition of connectedness: a set X C A is connected if any two sites in X can be connected by a path of links whose 
endpoints all lie in X. 

Returning to equation ( p2| ) , in order to recover the equivalent polymer system, we can rearrange it as 

Z(D)= fc M . ( 24 ) 

v : dv—0 

with (we miss the ininfluent factor N\) 

k(v)=e 2 ^^ J] 1/3 (v p ) . (25) 

pCA 

The main idea is to reexpress equations ( p4| ) and (|2^) in terms of closed 2-forms v with connected supports, which 
will become the supports of the polymers. With this purpose let us recall that it is possible to write [jl(| v — ^ Vi 
with the property that supp Vi are the connected components of supp v. Moreover one can write A v — 22 - A Vi with 
dA Vi — Vi and these relations allow the following factorisation: 

k(v) = H k(vi) . (26) 

i 

Finally, observing that with the above definition of connected set the condition dv — implies dvi — Vi, we can 
reorganise the sum in equation ( p4[ ) as follows: 

oc . n 

z ^) = E^t E tl K ( x i> D ) k(x,d)= Yl fc w • ( 2? ) 

n— Xi...X n i—l v: suppv—X 

The sum is extended to all finite connected subsets X{ C A, with the condition that X^ U Xj is disconnected if i ^ j. 
Thus we see that expression (^7|) has the form of a hard core interaction between polymers Xi of connected 
support. This allows us to use the techniques developed to deal with these systems in order to give the wanted 
bound on disorder field expectations. We point out that an expansion of the form given in (p7|) is common to many 
statistical (or quantum mechanical in the lattice formulation) systems: the polymer activities K(Xi,D) are the link 
to the original problem, depending on the form of the starting action. In particular the same expansion is possible 
for the U(l) gauge model with Villain action Q and gives monopole condensation in strong coupling fRM. The only 
difference from Wilson model is that the dual representation is gaussian and this simplifies the handling of polymer 
activities. 



B. Bounds on polymer activities 



It is a well known result that the main properties (mathematical and physical) of the polymer system can be taken 
in strict correspondence with the general behavior of the activities, which in turn depends on parameters such as the 
temperature or the coupling constant. In this subsection we give a bound on \K(X, D)\ which is known Jl3| to be a 
sufficient condition for the convergence of the Mayer expansion for log Z(D); moreover it will be of great importance 
in the estimate of [log Z(D) — logZ]. We want to show that VM>0 3 (3m such that for [3 < (3m 
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(i) \K(X,D)\ < e - M|x| 



'dD, 



■K(X,D)\<e~ M ^ 



(28) 



Here \X\ stands for the cardinality of the support of the X polymer. The basic tool used in the proof of equation 
(p8|) is the following upper bound on the modified (and rescaled) Bessel functions, 



/„(/?) <eW (£) 



'I 



(29) 



which easily follows from the power expansion Il5[ 







E 

k=a 



F(k + 1) T(n + k + 1) 



2A- 



The bound (29) tells us that for small /? the series whose n-th term is given by Ip{n) is convergent. Now let us sketch 
the argument that leads to (|2^). In what follows we assume that suppw = X: from equations ( p5| ) and ( p7|) one has 

(30) 



i#(x,£>)i< ^ n \h(v v )\<Y, n i^mi 

v : dv— pC^ a// v pCX 



We have added the contributions due to all integer valued 2-forms on X because we have in mind to exploit the 
convergence property of ^ I sin). In order to sum over all integer 2-forms with support on X, we first consider each 
set Yi of plaquette such that supply = X and sum over 2-forms v satisfying v p ^ Vp G Y- L . Then we collect the 
contributions due to all the sets Y{. Formally we have: 



\K(X,D)\< Yl Ell [m* 

Yi :supp Yi— X v.Vp^O p£Yi 



(31) 



Now let us focus on the generic term with fixed Yf. exchanging the sum with the product and using the parity 
property /„ (/3) = /_„ (/3), we obtain for the polymer activity 



\K(X,D)\< n 



Yi-.suppYi=X pGYi L n>0 



2 ]T I {n) 



(32) 



Thus to each plaquette in Yi is associated a factor that we extimate replacing Ip(n) by the upper bound ( |29| ) and 
summing the resulting series. Noticing that this is a geometric series missing the first term we are left with 



n e/3 < e -N p (Xi) log(3^) 



for j3 < - , 

e 



(33) 



where N p (Yj) is the number of plaquettes in Yj. From the fact that suppY, = X follows the relation Np(Yj) > j\X\ 
Moreover one can bound the number of Yj with support on X by e k ' x ' and obtains 

\K(X,D)\ < e -( A e- k )\ x \ 



* = \ l0S( ^ ) 



(34) 



Thus we see that part (i) of (|28|) is satisfied with Mp = Ap — k. As far as part (ii) of (28) is concerned, it can be 
obtained showing that there is a function G(f3) such that holds 



— 



K(X,D)\ < G(j3)\X\ 3 e-^- 1 ^ = F(\X\) ; 



(35) 



then one finds a constant Mp such that 



F(\X\) < e - M 'l x l . (36) 
For a proof of equation (p5|) see appendix A. The proof of (pq) can then be obtained by means of elementary analysis. 



G 



C. Bound on S2 (x, q ; y, —q) 



For sake of completeness let us now sketch the analysis given in Q to bound [log Z(D) — logZ]. Let us start from 
the definition of the Mayer expansion: 

00 ^ n 

logZ(D) = J2^ E MXi-Xn) 11 K(Xt,D) . (37) 

n=0 ' Xi...X n i=l 

Here ip c (Xi...X n ) is the connected part of the hard core interaction of the polymer system and it is nonvanishing 
only if Uj Xi is a connected set. The bound on the polymer activities allows the convergence of the Mayer expansion 
p3[ and this implies, for example, that correlation functions defined by differentiation of log Z(D) with respect to 
share the cluster property. 

Following the line of [10| we now define H(s) = log Z(sD) and observe that H(s) = H(—s) because the measure on 
equivalence classes [A is even. Hence we see that 



log Z{D) - log Z(0) = ( ds(l- s) H' 
Jo 



(s) , (38) 



where 

H"{s) = D ( b ) D (b') m(6, b') (39) 



b.iy 



<m') = E^E E 



MXi-x n ) x 



n=0 i,j = l bdB{Xi),b'dB{Xj) 



J] K(X kl sD) 



<k K{Xi ' aD) m K{Xi ' aDh m 



In (|40|) B(X) denotes the smallest rectangular parallelepiped in A which contains X. Moreover, by the exponential 
bounds on activities and derivatives given in equation (|28| ) follows that 



Yl K(X k ,.sD) 



Qp-KiXi, sD) —K(X h sD) < exp 



-Mj2\ X k\ 



(41) 



Since i/) c {Xi.,.X n ) =fi only if U, Xi is a connected set, we have that J2k \Xk\ > d(b, b') (d(b, b') is the distance between 
the two links) and so 



\m{b, b')\ < exp 



-Md(b, b') 



00 ^ n 

E ~\ E E 

n=0 ' i,j=l bcB{X i )y(ZB(X j ) 



ip c (Xt...X n ) exp 



\m e i* 



(42) 



Now, using the techniques of [|13j, one can work out from (42) the inequality 



\m(b,b')\ < exp 



-~Md(b,b') 



XDb 



-~M\X\ 



(43) 



where S is a constant little as fi decreases. For M large enough the sum over n appearing in (|43|) converges to a 
constant 8' which again is little as (3 decreases. Moreover reading the sum in (|3^) as scalar product between D and 
mD, we have the following bound on \H"(s)\: 

\H"(s)\ = \(D,mD)\ < \\D\\ 2 \\mD\\ 2 < \\D\\ 2 2 \\m\\ 2 . (44) 

From properties on matrix norms we have that ||m||2 < supv (^2 b \m(b, b')\): using the previous result we obtain 



\H"is)\<S'\\D\ 



<P{P) 



(45) 
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with p((3) — > as (3 — > 0. It is important to point out that the argument works because ||-D||2 is bounded uniformely 
in x and y. Finally, using (^) and (^9|) we conclude 



S2(x,q; y-q) > exp 



1 



(46) 



This relation implies in particular that Si(x, q), defined by the limiting procedure in which y — > oo, is nonvanishing. 
In the language of field theory we can say that the field describing magnetic monopoles aquires a nonvanishing vacuum 
expectation value. This implies the spontaneous breakdown of the topological symmetry associated to the magnetic 
charge conservation and signals confinement of electric charge [Q . 



III. MONOPOLE SECTORS IN THE QED PHASE 

In this section we prove the relation (|l^) for the soliton two point function, showing that for (3 large enough and 
\x — y\ — > oo it is possible to find a positive constant m{(3) such that 

S 2 (x,q;y,-q)<e- m ^ x -v\ . (47) 

In order to proof our statement, we use a slight modification of the expansion given in jll],[l2j and p6| : we reexpress 
the partition function as a gas of monopole loops, to which apply a renormalization transformation. Estimates on the 
renormalized loop activities enable us to extract the relevant contribution to £2 (a;, q ; y, — <?)■ 



A. Expansion of Z(D, Q.) in interacting monopole loops 



From equation (^), defining the modified partition function in dual representation, it is natural to introduce a 
measure on equivalence classes [^4] of integer valued 1-forms given by 



d»(A) = ± Yl I {dA p ) f . 

Z p ca J W 



d/i (A) = 1 



(48) 



The main idea |Tl]] on which is based our construction is the following: we want to introduce a measure dfj,i ff (A) on 
R n (n is the number of links in A) , which reproduces (Q) once we constrain the real variables A^ to integer values 
and pick them on a gauge slice. Such a measure should enable us to make suitable extimates in the weak coupling 
region. We can fulfill the first constraint inserting a sum of 5- functions for each link variable A^, which by Poisson 
summation gives us the monopole currents; formally we have 



dli{A) = \ n 



6CA 



E 

q&Z 



; i27T q A b 



d[ii„ (A) 



(49) 



Moreover, since we are going to compute expectations values of gauge invariant observables, the only contributions 
come from conserved currents; hence we can impose 

(50) 



(51) 



dni/, (A) e l {J ' A) =0 for S J ^ . 

Actually it is possible to construct a measure with the required properties, taking the limit A — > Z of 

dm, {A) = J] T P ^ H dA b A b eR . 



beA 



1/3 {4>) is a suitable interpolation of the modofied Bessel functions, which has been constructed in (T^| for large values 
of P, and has the properties listed below: 



1. 



1 f 27r 



d9 e 0cos6 e l4 



for integer 
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2. ((f>) is an even, positive and integrable function on R; 

3. I/3 (4>) is analitic on the strip |Zm</>| < ^. Moreover it is possible to find a constant c such that 



,1 



P 1 



■ + ia) . 



< exp 



with < 5(a) < ce 



27r|a| 



for 1< \a\ < - . 
- I I - 2 



With these new tools the soliton two point function can be written 



(52) 



Henceforth, for notational convenience, the factor 2?rq will be absorbed in D and f2. Our first objective is a suitable 
rearrangement of the product appearing in (|52"|), which, by the constraint (|5^), gives the usual coupling of A with 
external current loops. 

Let us start defining a current density p as a 1-form on A with values in 2nZ. An 1-ensemble £ is a set of current 
densities {p} whose supports are disjoint and such that dist(p,p') > 2? Vp, p' £ £. It is useful to collect the above 
currents in 1-ensembles using the following property p2[ : 



n 

6eA 



i2ir n b A b 



E 



n 



+ cos(A,p)] 



(53) 



The index a runs over a finite set, each £ a is an 1-ensemble and d a > 0. Moreover, the bare loop activities satisfy 
< K(p) < e' "''" 1 , where \\p\\i is the norm given by ||p||i = ^2 bf z suppp \pb\- The next step is to consider the string 
Q as a current density and construct suitable 1-ensembles containing open-ended currents which are obtained by 
grouping the currents "touching" f2, in a sense which we'll specify below. We choose f2 (by hyper-gauge invariance) 
such that suppfi R suppD = 0: This means that we are taking f2 with support in the region bounded by the 
hyperplanes z a = x° and z° = y°. The disorder field expectation can be written as follows |ll| , |l2] , ^6| : 



Z(Cl, D)= J dp If) (A) c T [cos{A, D~n)+ K(p T ) cos(A, D - Q. + p T )} x [1 + K{p) cos(A, p)] 



(54) 



Now, some comments on equation ([54]): the currents p T are divergenceless and their support has non vanishing 
intersection with N = {b £ A : dist (b, supp (D — SI)) < 1}: p T are the currents touching Q. c T are positive constants 
and K(p T ) satisfies < K(p T ) < e' c " pT " 1 . Moreover £ T is an 1-ensemble of divergenceless currents having vanishing 
intersection with N. We must point out that £ T U {p T — 51} is still an 1-ensemble: we have divided the closed currents 
from the open-ended p° T = p T — ft, which play a peculiar role in the proof of clustering of the soliton correlation 
function, as we shall see below. 



B. Renormalization transformation 

Now we make on ( |54| ) a transformation that renormalizes the activities of currents and allows us to give the wanted 
bound on S2(x, q ; y, —q). The transformation consists in the explicit integration of the factors containing the currents 
p on a suitable subset of suppp, which we call B p , characterised by the fact that two different links contained in it 
belong to different plaquettes and 

E \Pb\>^\\ph ■ (55) 

In dimension four the geometric constant c can be fixed to be ■ Such a renormalization is based on the application 
of the following property. Let us consider a function G( A) which does not depend on .4;, for some link b £ A: for 
arbitrary real a such that \a\ < 4 then we have 

G(A) dm, (A) , (56) 



e tqA " G(A) dm, (A) = e - M ^ a ' g) / e lqA 



W i l 3{e{b lP )a,dA) 



p:b£dp 
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where 



Tib 



Ep(a,q) = qa- — g(a) and ip(a,4>) 



+ ia) g(°) 
e "T - 



(57) 



In equation ( |57| ) rib is the number of plaquettes containing b and e(b,p) is a factor ±1 which gives the orientation of 
b in dp. The proof of this property is immediate if we make a complex translation on the variable A ( A — > A + ia) 
and use the properties of the function Ip (</>) . 

Now let us focus on the generic term in the sum over the index r appearing in ([si]). First, we choose S p C supp p for 
every p G £ T and £> pT c supp p^ . By the defining property of these subsets and the fact that we are dealing with 
1-ensembles, all links selected in this way belong to different plaquettes. After exponentiating the cosines in (p3) and 
decomposing the resulting products, this observation allows us to apply the relation ( p6| ) to the resulting terms for 
all links in B p , B p -r and suppfi. These integrations produce the transformation of cos(yT, p) into a function cp(A, p) 
and renormalize the activities K(fl) as described by the following relations: 



Z(fi,D)= ( dp h (A)J2^ [z(ft)cp(A,D-n) + z(ft,p T )cp(A,D + p n -)] x J] [1 + z(/3, p) C/3 (A, p)] 
The renormalized activities are given by 



(58) 



n 



supp Q, 

z(ft,p T )=K(p T ) Yl e -S (a,(D+p^ )h ) . 



beer 



z(p )P )=K(j>) J] e -^ (a ' pb) . 
The renormalized version of the cosine is (A, p) = 72.e [e^j (A, p)] where 



(59) 



JJ ^ (ea, cL4(p)) 

p6T(£? p ) 



and 



e/3 (A,p) = e^) 
From the relations given above it is clear that \cp {A, p)\ < 1. 



T(B P ) = {peA : be dp b e B p } 



(60) 



C. Estimates on renormalized activities and cluster property 



In order to extract a bound that assures the cluster property we must now choose a suitable value of the parameter 
a appearing in Ep. By property 3. of the function Ip (</>) follows that 



5/3 (a,q) < e -Ep (a,q) 



with E (a,q) = qa- ^ce 2l|n| 



(61) 



In order to give an upper bound on activities as strong as possible, we take the value af a maximizing Ep(a,q) in 
the domain \a\ < £ and we denote it by a rn . For a fixed value of ft it turns out that a m depends on the value of the 
parameter q, which stands here for the value of the generic current on a link in B p (q — > p& in equations (]5^)). We 
find 



and correspondingly 



2lT \ l-KUbC J 



if \q\ < m 



°>m{<i) = £ if \q\ > fa 



ft 



E™ (q) = Ep(a$,q) 



m 

2n 



log 



( m 

\2irribC 



- 1 



Ef ( q )= E(3 (a£\q) = \q\^ 



_ 2ntc p 
\q\ft 2 



(62) 



(63) 



The discriminant value q is defined by a m (q) = The important feature of these current self-energies exctracted by 

renormalization is that for ft large enough both E^ (q) and E^ (q) are positive. 

With the above choice for the parameter a the renormalized activity of the generic current p satisfies 
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,(p, P )<K{p) n 



-\ Pb \A(0,p b ) . 



beBr, 



1 

l«0 



log f #M) 



2n b c tt/3 



WW 



for |p b | < g 
for \p b \ >q 



(64) 



In both cases A(f3,pb) is bounded from below by a function which does not depend on p\>: 



A(f3, Pb ) > 



1 

2tt 



logf^ 

iU & ^ 127TC 



- 1 



= (/?) 



I i - 40) = a^ w 



(65) 



The first bound is obtained using the inequalities \pb\ > 1 and n b < 6; the second one is obtained replacing \pb\ with 
q. We must piont out that both A^ (f3) and A^ 1 (/3) are positive functions increasing with (3. If now we define 
A(P) = minjA^ 1 ^ (P),AW (0)}, using the properties of B p and K(p) we can write 



z(p, p) < K{p) e- £A ^ II"" 1 < e"(^W- fc ) Uz-lli 
In particular one can bound the renormalized activity of the string as follows: 

z(j3) < e" 5A(/3) H " 1 < e -ZM0)i\*-y\ . 



(66) 



(67) 



More delicate is the estimate of z(f3,p T ), because of the presence of the Coulomb field D. In fact one has B p t C 
supp(p^ T ) but the generalized current density is (p Qr + D). In general for these activities holds 



The starting point to prove (p8h is the relation 



,(P,pT)<K(pT) Yl e 
beB„T 



-MP)\pl >T +D b \ 



(68) 



(69) 



First one can easily see that in the case in which supp p T R supp D — one can deal with the current p Q ' T = p T — il 
as for the common p (see equation (66)) and from the relation Ip^H >q\x — y\ follows d68| ) with C(f3) — 1. 
The case in which supp p T PI supp D ^ can be worked out with the following trick: we distinguish the links b 6 B p r 
such that \D b \ > ^ from those such that \D b \ < i (in other words we decompose B p -r = B^-r U B pT ) obtaining a 
factorization in equation (|69|). Noticing then that p T takes values in 2irZ, for b S B% we have 



\pb T +Db\ = \pf\ 



D 



l + ^)l>l^li 



9 h 



(70) 



by which we reduce the factor involving B pT to the standard form (|64|). The term involving B p T will give us the 



constant C(/3). In fact after little simple algebra we obtain 



z(P,p T ) < K{p T )e-i £AW ^ SlT ^ G{D,p n -) 



with G{D,p n ^) = Yl exp 



Am 



\pf 



\pb 



Finally it is possible to show that G(D,p T ) is bounded by a constant C((3) only dependent on (3 and the number of 
links on which \D\ > i (actually this number is a function of q); this completes the proof of equation (|6S|) . 

Now we have all elements to complete our proof. By the property that \zcp\ < \z\, starting from (pq) we can write 



\S 2 (x, q ; y, -q)\ < sup T {z{0) + z((3, p T )} x / dp I(i (A) ^ c T J[ [1 + z(0, p) c p {A, p)] 



(71) 



Let us call TZ the integral in (fry): considering the explicit form of the normalization factor Z, it can be represented as 



K 



Er Ar 



(72) 
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with 



A T = c T f dpi, (A) H [1 + z(J3, P ) c fj (A, p)} ; (73) 

B T = c T f dp h (A) [1 + z[fi, p T ) c (p\A)} 11 [1 + z(fl, p) c (A, p)] . (74) 

It is simple to see that for j3 and \x — y\ large enough we have A T < 2 B T . In fact choosing suitable values of (3 and 
\x — y\ we can obtain that z(j3 1 p T ) < \ and z(/3, p) < 1; thus in order to evaluate B T we must integrate in the positive 



measure 



dpi, {A) II [l + KAp)c/)(Ap)] 

the function [1 + z((3, p) cp (A, p)] > i. This observation leads to the conclusion that 1Z < 2. Finally using this result 
in (71) we get the relation 

\S 2 (x,q;y,-q)\ < 2 sup T {z{(3) + z{(3, p T )} . (75) 

Cluster property of the soliton two point function follows then by estimates of equation (^) and (|68| ) on renormalized 
activities. 

IV. CONCLUDING REMARKS 

In conclusion we summarize the results: in the weak coupling region the Hilbert space of the reconstructed Lattice 
Quantum Field Theory splits into orthogonal sectors labelled by magnetic charge. Instead, in the strong coupling 
region the lattice solitons do condense in the vacuum sector: the symmetry associated to the topological conservation 
of magnetic charge is spontaneously broken, as signaled by the nonvanishing expectation value of the charged monopole 
operator. Indeed this is just the criterion for quark confinement proposed by 't Hooft M. 



Moreover our analytical results are in agreement with numerical calculations |17|: these show that the parameter 
regions in which S\ > and S\ — coincide with the confining phase and QED phase, respectively, and allow to 
extract information about the behavior of the system at the transition. Our result suggests that the correlation 
functions of soliton (disorder) fileds are indeed a viable tool for the study, both numerical and analytical, of phase 
transitions in lattice models that exhibit monopole-like topological excitations (we point out that a slight generalization 
allows the extension to lattices with non trivial topology) . Along this line there is the possibility to study analitically 
the disorder fields associated with vortices in the 3d XY model and, although more remote, a generalization to non 
abelian models. 
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APPENDIX A: BOUND ON ACTIVITY'S DERIVATIVES 



In this appendix we sketch the argument which leads to equation (|35|) . The proof is based on the following inequality, 
given in (lOl: 

\\A v \\<^-{v,v) 2 <cN 2 p m A v , (Al) 
where N p is the number of plaquettes in X and 

m v = maxp C x \v p \ . 
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The derivative with respect to D b modifies the expression of k{v) by the moltiplicative factor 2tt A v (b) and thus we 
have: 

' } K{X,D)\<cNlYl k * E UhM (A2) 



dD b 

k>0 [v]:m v =k PCX 



Now we can write 



E II^(«p)= E E U^m- ( A3 ) 

[v]:m v =k pCX Yi:suppYi=X [v]:m v =k,v p ^0 peYi 

The 2-forms with m v — k may be obtained fixing the value of v to ±fe on a plaquette pm and summing over 

a k 

configurations of integers n < k in the remaining (Np(Yi) — 1) plaquette. We extract a factor ^ (coming from 

Ip(v PM = k)) from each term with fixed pu and notice that the contribution due to all configurations is bounded by 

e -(N P (Y z )-i)A^ The resmt is 



" K(X,D)\ <cNl fc 4 e-W- 1 )^ Yl 



fc>0 v ' Yi:suppYi=X 



(A4) 



The series in k can be estimated by the value of the integral of x 4 e Bx on R + (B = log(^)), and this gives G((3) 
apart from factors. Moreover the sum over Yi is bounded by 2|X|e fc ' x ' in such a way that we can write equation (|3 



\JLk(X,D)\ <G((3)\X\ 3 e-U x \-V M *=F(\X\) . (A5) 
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